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APPROXIMATE ANALYTIC SOLUTIONS FOR THE OPTICAL 
PUMPING OF FLUORESCENT DYES 


Nabil M. Lawandy 
Goddard Space Flight Center 
Greenbelt, Maryland 


IntrcKluction 


I>ye lasers offer spectroscopists tunable and coherent radiation from 33<> 
to 117<> nanometers (nin). Tunability is achieved by manipulating dye concen- 
trations, solvents, and resonator losses. Laser-pumped pulsed dye lasers are 
especially useful for time-resolved spectroscopy Ixjcause the laser pulses have 
the same shape as the pump pulses above threshold. However, pulses pro - 
duced by those lasers arc not always ;is energetic as might be required. In 
general, the stringent requirements on spectral purity limit the pulse energies 
to a few hundred microjoules. 

A possible remedy to the low pulse energies of such systems is regenerative 
or multipass amplification using flashlamp-exited amplifiers [1], [2], and [3], 
The theoretical imtler stand big of such :ui amplifier is based on a knowledge of 
the time -dependence of the amplifier inversion. Because of the order of mag- 
nitude variance in flashlamp rise times and triplet-state deactivation rates, 
triplet-state accumulation must lie treated as a variable. Sorokin et al. [4] have 
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produced computer solutions to tlu* rate equations for short-pulse laser pumping. 
Schafer and Sehinidt | a| have given analog eomputer plots for the excited singlet- 
state population in the limit of very long triplet lifetimes. Moreover, Keller 
|*i) has examiniM the effects of triplet states on the efficiency of a dye laser using 
computer solutions and neglecting the depletion of the ground-state population 
due to pumping, lie also provides steady-state solutions for the nuntbor of 
photons in the laser cavity. 

1'his paper develops a technique for solving coupled-rate equations to a high 
degree of accuracy. The method is applied to the inversion in a flashlamp- 
pumped dye solution. A three-level approximation is used to describe the dye- 
field interaction. 

Dyc-I* umpmg Mechanisms 

A large molecule, such as a typical dye, has many normal vibrations that 
arc coupled to electronic suites. Transitions between such states are broadened 
by collisional and electrostatic perturbations with the solvent molecules. More- 
over, each vibronie sublevel of an electronic state lias rotational fine-structure; 
however, this is smeared out by the rapid rate of collisions (10 ** collisions/s). 

The pump radiation in this system corresponds to an electronic transition 
from the ground vibronie level to an excited vibronie level. This leads to a 
nonequilibrium Frank-t'ondon state. Furthermore, the transition is between 


singlet states and is strongly allowed (1 in figure 1) The state, Sj, quickly 
thermalizes (10 -1 “ s) because of solvent collisions. At this stage, a molecule 
may either fluoresce to the lowest electronic state (2 in figure 1) by Frank- 
Contlon selection rules or make a nonradiative "intersystem crossing" to a lower 
triplet state at a rate, k^-p. This is detrimental to laser action in three ways: 

(1) it competes with fluorescence; (2) it creates optical losses by placing mole- 
cules in triplet level Tj, which can make triplet-triplet transitions (Tj—Tg ln 
figure 1); :uid (3) it depletes the ground-state population and therefore the 
number cf molecules in the fluorescence cycle because of the long lifetime of 
Tj. Phosphorescence from T. to the ground vibronic manifold is spin-forbidden 
and may trap molecules for limes ranging from milliseconds to seconds. How- 
ever, the use of paramagnetic collision partners in the solution (Oo for example) 
that induce a strong spin-orbit coupling can decrease the triplet lifetime through 
radiationless deactivation. In addition to the previously mentioned processes, 
singlet-singlet absorption in figure 1) may take place and cause additional 

optical losses. These are not as significant as triplet-triplet losses because of 
the relatively small population of Sj. 

Using this description for the interaction, the system can be modeled with 
rate equations and phenomenological decay rates. If the thermalization rates 
are approximated within each singlet manifold as being instantaneous, the system 
may be described by three levels: S 0 , Si, and Tp In the absence of stimulated 
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Figure 1. Energy level representation of a dye molecule. Lines 1 and 2 denote 
Excitation and fluorescence, respectively. Dashed lines indicate non 
radiative transitions, and solid lines indicate radiative transitions. 
Intersystem crossing is described by a rate, k^., and the lifetime 
of the triplet state by T . 




emission, the populations of the levels obey the following diffcr*ntial equations: 
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where l(t, \» is the flashlamp intensity as a function of time, and g(\> is some 
line-shape function of the- pump wavelength, \ , 
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Equation 5 is a matrix equation with the matri/.ant solution given by: 
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whore N(0) 

0 

Is the initial population vector. 


V 

• 


0 



r> 


T 


because M(t) is not in diagonal form, ve can expsind the exponential matrix 


operator as a power series. This yields: 


KXP^ / M <t>dt] I ♦ J M(t j)iitj f 1/2 j ( M (t 0 )M(t.j)dt 1> dt, J • ... 


For tlie particular matrix of concern. 

-yt 0 Z(t) 

fittft = k ST t -t/T T 0 

t/T t/T -Z(t) 
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where 

Z(t) /\V(t)dt (j>) 

This formulation can be used to approximate the ground-state population, 
N () (U, to an arbitrary degree of accuracy. The ground state is chosen for two 
reasons: (1) it is the least sensitive to the magnitude of its own first derivative, 
and (2) it is a function of both N ^ (t> and N, ( .(t). The approach to solving for X^t) 
is based on solving the first-order differential equation (la) using an approxi- 
mation for X ( (t) derived from the expansion given by equation 7. Designating 
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the order of the N () (t) approximation hy a superscript, 

V» • N o 


N o (t) = M oP ' Z(t, J 

(ID 

2 N 0 (t) - N Q [1 -Z(t)(l - t/T s - Z(t>] 

(12) 

Similarly, for N,j.(t) 

°N. r (t) 1 N T (t) = 0 

(13) 

V r (t) - N 0 k ST ltZ(t) l 

(14) 


The first-order expression for N () (t) can be used to solve equation la for 
Xj(t). Examining equation la, we see that it has the standard solution: 


Nj ft) 
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EXP(yt) N <t)W<t)dt 



(15) 


where C is chosen so that N (0) 0. Substituting equation 11 into equation 15 

gives: 


Nj (t) = 


N n EXP(-yt) 


/iSXP(yt) ^1 - Z(t)J W(t)dt + C 


(1(5) 


This expression accounts for the depletion of the population hut assumes that 
all molecules entering T ^ remain there. 
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Finite Triplet lifetime 


The accumulation of molecules in T combined with a relatively small value 
for T,^ can rosult in a sizable return rate to the ground level. Using tlu; expan- 
sion (equation 7), we arrived at a value lor N^,(t) given in ociuation Li. The 
triplet contribution to N () (t) appears in equation lc as Substituting 

equation 11, the contribution to N () (t) by this rate can be estimated by the term: 


Nk r 

0 ST / t/(t)dt 

T t J 


Equation 17 e:m lx‘ used to modify either equation 11 or equation 12 before using 
them to solve equation la. The modified term is given by? 


1 N 0 (t) S Q |^1 - Z(t) ♦ k gT /T T j tZ(t)dt 


Linearly Hi sing Flnshlanu 


This section examines the approximations derived for the ease of a pump 
source with a constant slope, which has the form: 


l(t,A ) mt 


Equation 1!) results in: 


Z(t) Il(A)mt /2 
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Substituting equation 20 into equation 10 amt integrating yields for the population 


of S ? 

N ! (t) N 0 ll(\>m/y t 

where 


l/y - B(x)mV(t)/2 * EXP(-yt) 


l/y - 3B(x)m/y 



V(t> (t‘ l - 3t“/y ♦ utA 2 - G/y' 1 ) 


( 22 ) 


Expression 22 can he modified to account for molecules that make the transition 
from Tj to S >( l>y adding the term? 


N 0 B(\) 2 m 2 k ST R(t)/T T r 


(23) 


where 

K(t) 


[t ’ - 5t l /y * 20t'Vy“ - (i0t 2 /v 3 ' 120t/y 1 - 120/y’* ► EXP(-yt)/y J ] (24) 


Expression 21 modified by equation 23 is plotted for various parameters in 
figures 2 through I. Figure 2 shows the effect of varying the pump slope. Be- 
cause a relatively long T was chosen, flashlamps that rise quickly produce 
higher peak populations but are short-lived since most of the molecules become 
trapped in T . However, for smaller pump slopes, the triplet return rate is in 
Ijctter competition with pumping. This permits the population to bo maintained 
for longer periods at lower levels. 


Figures 3 and 1 show that k and T play complementary roles. The 

i) 1 1 

duration of inversion is strongly dependent on l>oth rates. As figure 4 indicates, 
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* ^ ure 3 * l ,( »Pulation density of Sj as a function of time for two values of the 
intersystem crossing rate k sr . The relevant constants are. 
m 3 * lO^l, 5 „ 10 -9 8 , and T ^ x x 1( T fi 8 . 
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figure* I. Population density of S as a function of time* for two values of the 

* 1 1 
triplet lifetime, I' j . i’he relevant constants are? m It * 10 , 

-9 . - 7-1 

I ft * 10 s, and k r* o 10 s . 

s ST 

■ 
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for T values that are sliort (~k >, we can achieve a nondumped inversion. 

I c) 1 

Depending on the resonator losses, such a situation could result in continuous- 
wave laser operation. In actual practice, the shortening of T,^ is accomplished 
by using some quenching agent. However, most quenching agents are not 
specific for triplet states and therefore also increase k . On the basis of how 

u 1 

this affects each rate, a net enhancement may or may not Ik* achieved. 

The analytic approximations given exhibit all jie coarse features of the 
population of as a function of time. The plots given agree well with respect 
to both peak population density and time-dependence when compared to computer 
solutions. For comparison of the two solutions, the reader should consult 
deferences 4 through (i. The results given can be readily used to estimate the 
single-pass gain as a function of time for a particular dye. Moreover, the same 
technique can be used in conjunction with a fourth rate equation for the cavity- 
photons to examine transient laser action. 

Conclusions 


This paper describes how the use of a matrix expansion can lead to analytic 
approximate solutions for a component of a set of rate equations. The application 
of this approach to the optical pumping of dyes illustrates the validity of the* re- 
sulting expressions. All the important features concerning the dye problem can 
be found in the first-order solutions given. More subtle effects leading to 
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ik' milnl structure ran he found by including higher order terms in the matrix 


expansion. 
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